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Abstract 



The possibility of a type of semiconductor quantum dots obtained by spa- 
tially modulating the spin-orbit coupling intensity in III-V heterostructures 
is discussed. Using the effective mass model we predict confined one-electron 
states having peculiar spin properties. Furthermore, from mean field calcu- 
lations (local-spin-density and Hartree-Fock) we find that even two electrons 
could form a bound state in these dots. 

PACS 73.21.La, 73.21.-b 



1 



Typeset using REVT^ 



Low-dimensional semiconductor structures constitute one of the most suitable scenarios 
to observe quantum phenomena in condensed matter physics. In this sense, the atomic-like 
properties of quantum dots and the conductance quantization in narrow wires are semi- 
nal examples widely studied [1,2]. At present, the available fabrication techniques of such 
semiconductor structures rely on charge-based confinement at length scales where quantum 
effects arc unavoidable. More precisely, confinement is obtained through the energy band 
discontinuities of different materials or by directly creating a potential landscape by means 
of metallic electrodes on a semiconductor heterostructure. 

The aim of this work is to point out, using a theoretical model, the feasibility of a class of 
mesoscopic semiconductor structures where the spin plays a relevant role in the mechanism 
that determines the confining properties, without being a pure spin-dependent confinement 
[3]. The mechanism also differs from conventional electrostatic confinement in that it does 
not rely on the direct coupling of the electron charge with a horizontal electric field. 

In III-V non-magnetic semiconductor heterostructures there are two relevant interactions 
concerning the spin: the spin-orbit coupling and the hyperfine interaction between carriers 
and polarized nuclei. The latter mechanism has a characteristic energy scale usually lower 
than that of the spin-orbit coupling [3]. 

In semiconductors, spin-orbit coupling stems from the relativistic effect caused by the 
electric field due to the lack of inversion symmetry of certain alloys like III-V heterostruc- 
tures. Depending on the particular origin of the electric field two contributions are distin- 
guished; the electric field created by the bulk inversion asymmetry (BIA) of the material 
gives rise to the Dresselhaus term [4], while the asymmetry in the profile of the heterostruc- 
ture, structural asymmetry (SIA), generates the so-called Bychkov-Rashba term [5]. This 
latter mechanism constitutes the basis of the proposed electronic nanostructures. 

As mentioned above, the intensity of the Bychkov-Rashba spin-orbit coupling depends 
on the effective electric field perpendicular to the plane of the quantum well. This implies 
that it can be, in principle, engineered by adjusting the specific fabrication settings of 
each heterostructure. Another possibility to control the Rashba interaction, that has been 
demonstrated experimentally [6-8], is by means of an external electric field perpendicular 
to the plane of the quantum well. 

Exploiting the tunability of the Bychkov-Rashba mechanism we consider quantum dots 
whose main ingredient is a space-modulated spin-orbit intensity. This spatial modulation 
could be obtained with biased small electrodes on top of the semiconductor heterostructure 
[9,10]. However, such an arrangement would not only modify the vertical electric field but 
it would also generate an in-plane gradient of electric potential, leading to hybrid quantum 
dots where both mechanisms are present: horizontal electric potential gradient and spin-orbit 
intensity gradient. A more refined method to obtain space modulated spin-orbit intensities 
requires a deeper insight in the origin of the Rashba intensity. In this respect, de Andrada 
et al [11] discussed the relevance of the barrier layers in the determination of the spin-orbit 
intensity. More recently, Grundler [8] has shown that for an InAs-based quantum well an 
applied vertical electric field of the order of 10^ V/cm can lead to a noticeable control of 
the spin-orbit intensity, whereas the built-in electric field of the sample is of the order of 
10^ V/cm. This high sensitivity to a relatively weak electric field has been attributed to the 
induced variations in the wavefunction penetration into the barriers. The strong relevance of 
the barrier layers in the determination of the Rashba intensity supposes a great advantage 
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that could help to create a quantum well having a space-modulated spin-orbit coupling 
through space-engineered wavefunction barrier penetration (using inhomegeneous barrier 
layer composition, thickness, ...). This would lead to a fixed (not tunable) substantial step 
in the spin-orbit intensity capable of inducing the confinement. 

We model structures confining electrons from the conduction band of the semiconductor 
and we use the effective mass approximation. We also suppose that the motion perpendic- 
ular to the well's plane is restricted to the first transversal subband, so that the electrons 
effectively move in a two-dimensional region. These simplifications might be quite drastic in 
some cases and, therefore, the scope of this work is hmited to a discussion of the feasibihty 
of the effects rather than to a detailed quantitative description. In this model we include 
the standard Bychkov-Rashba term 



h 



(1) 



where Px and Py represent the components of the in-plane electron's momentum and the cr's 
are the corresponding Pauli matrices. 

In Eq. (1) the inhomogeneity of the spin-orbit intensity is included through the spatial 
dependence of the coupling constant Aij(r), which is assumed of radial type (r and 9 are 
used to label the standard polar coordinates in the plane) . To characterize the quantum dot 
we consider the following step-like variation of the spin-orbit intensity 

Mr) = Ae + (A, - Ae)- ^ , (2) 



1 + e- 



where Ae and A, represent the external and internal constant values of spin-orbit intensity, 
respectively. The quantum dot radius is given by Rq and cr is a small diffusivity introduced 
to avoid discontinuities in the numerical solutions. 

The Bychkov-Rashba interaction breaks the rotational symmetry of the Hamiltonian. 
Nevertheless, a symmetry combining spatial and spin degrees of freedom is maintained when 
the Rashba intensity is a purely radial function. Namely, the Hamiltonian is invariant under 
a combined rotation in real and spin spaces by the same angle, i.e., it commutes with the 
generalized angular momentum = + S^- Note that, at variance with conventional 
circular quantum dots (based on electric-potential confinement [1]), the Hamiltonian no 
longer commutes with and Sz separately. As a consequence of the J^-symmetry the 
single-particle solutions can be expressed, in a general form, as two-component spinors of 
type 



0„,i(r)e^O-+V2)^ 



(3) 



where n and j are quantum numbers characterizing the states, with values n — 1,2, . . ., 

j — ±1/2, ±3/2, Since the angular parts are analytical, we can reduce the Schrodinger 

equation to one-dimensional radial equations for the two spinor components. 
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where 0^ and (f)i denote 0niT(^) ^-^^ 0njj,(?'), respectively. 

The above model has been analyzed using two alternative numerical procedures: a) 
solving the coupled one-dimensional equations for the radial components of the eigenspinors, 
Eqs. (4); and b) solving the full two-dimensional problem in Cartesian coordinates, without 
imposing any symmetry restriction, in a uniform grid. Before discussing the numerical 
results it has to be noted that systems with the above spin-orbit interaction are still invariant 
under time reversal. Consequently, a two-fold degeneracy, known as Kramers degeneracy, 
must hold. From the radial equations it can easily be seen that the substitutions: 0^ 4'ii 
4>i ~^ a-iid j ~^ ~j Isad to a new state (Kramers conjugate) with opposite generalized 
angular momentum and the same energy eigenvalue {Snj — Sn-j)- The expectation values 
(Sz) and (Lz) also change sign with the Kramers transformation. It is worth stressing that 
in the usual quantum dots (with an electric potential confinement of circular symmetry) the 
degeneracies are 4 and 2; nonzero (zero) angular momentum states being four-fold (two-fold) 
degenerate. Therefore, the reduction in degeneracy from 4 to at most 2 is a genuine effect 
of the spin-orbit interaction. 

Since the spin orientation in Eq. (3) depends on the position r, the eigenstates will show 
characteristic spin textures. The distribution of vertical and parallel spin read 

(a.)(r) = |0.,t(r)|2-|0„,i(r)|2 
(a||)(r, e) = 2(l>r,j^{r)(j)nj^{r) [cos{e)x + sm{e)y] . (5) 

Note that while the vertical component is a pure radial function, the in-plane (horizontal) 
spin points in the radial direction. The fact that the horizontal spin follows the position 
vector can be easily explained from the combined rotational invariance in spin and position 
space: consider two points in the dot having the same radius and different azimuthal angle 
9] due to the symmetry of the Hamiltonian with J^, the local in-plane spin orientation in 
the second point must be that of the first point rotated the angle needed to go from the first 
point to the second in the real space. This property of the eigenstates is of particular interest 
for determining the role of the geometrical phase (Berry's phase) in mesoscopic structures 
[12,13]. 

Prom the numerical simulations we observe that a necessary requirement for confined 
states is that the spin-orbit intensity within the quantum dot must be larger than that of 
the outer region, i.e., Aj > Ag. The reason for this behaviour can be understood from bulk 
considerations: the bulk energy bands 

display a constant offset —m*Xl/2h'^ that depends on the spin-orbit intensity (for a schematic 
representation of the bulk bands see Fig. 1). Therefore a bulk uniform sample with a big 
value of spin-orbit intensity has an energy origin for the eigenstates lower than that of a 
sample having a smaller intensity. In a finite system this energy mismatch can lead to bound 
states in the region of higher intensity. Nevertheless, it has to be noted that for each spin 
orientation the energy mismatch is only confining for a certain range of momenta. This 
means that the confining mechanism is spin-selective through the motion of the particles. 
In the case of circular symmetry, it implies that electrons rotating in one sense are confined 
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if they have a given spin orientation, while electrons rotating in the opposite sense should 
have the opposite spin orientation to be confined. 

We stress that the SO-dependent ofi^set in Eq. (6) is the real cause of confinement since 
this term acts as an electrostatic potential when the SO intensity is varying in space. Ac- 
tually, a different SO Hamiltonian in which this offset were not included would not lead to 
confined states. 

Another question of relevance refers to the magnitude of the step in the spin-orbit inten- 
sity needed to obtain confined states. For a typical size of dot radius of the order of 10 atomic 
effective units [14] simulation shows that an intensity in the range Aj ~ 1.5 — 2.0 x 10~^ 
eVcm is required to obtain a single-particle bound state when the external intensity is zero 
(Ae = 0). It is also important to specify the external value Ae, since the minimum step 
required to achieve confined states strongly depends on this quantity. For instance, when 
the external intensity takes a value of the order of Ae — 2 x 10^^ eVcm the step needed is 
reduced and the required internal value is in the range Aj — Ag ^ 0.4 — 0.7 x 10~^ eVcm. 

Quantum dots satisfying the required conditions stated above show a level structure 
analogous to that of conventional quantum dots. For instance, a dot characterized by spin- 
orbit parameters in the range of the biggest Rashba intensities experimentally measured 
[8], Aj ~ 4 X 10~^ eVcm and Ae ^ 2 x 10~^ eVcm, shows a total of 14 single-particle 
bound states distributed in two subbands (see Fig. 1). The first one containing 10 states 
characterized by the quantum numbers n— 1, j = — 9/2,— 7/2,...,7/2,9/2 and the second 
one having the remaining 4 states n — 2, j = —3/2,-1/2,1/2,3/2. Figure 2 shows the 
spatial representation of the lowest energy state of this dot corresponding ton = l,j = 1/2 
and its conjugated state n = l,j = —1/2. The upper panels show the radial profile of the 
spinorial components for both states, the only differences between the two cases are the 
exchange of the spin index and the sign change in one of the components, in agreement 
with symmetry properties derived above. Middle panel shows the spatial distribution of the 
density which has radial symmetry and is common for both conjugated states. Finally, the 
lower panel displays the radial distribution of the 5*^ spin density for the two states, showing 
that conjugated states posses opposite spin characters. 

Up to now, we have only presented results concerning the single-particle properties of the 
dots. There is also the possibility of the spin-orbit-based confinement to be robust against 
electron-electron interaction and, thus, able to support multielectron bound states. To treat 
the interacting problem we impose no symmetry restriction in space and we solve numerically 
the full two-dimensional problem including the interaction within density functional theory, 
using a generalization of the local-spin-density approximation (LSDA) for non-coUinear spin 
densities (see Ref. [15] for a previous application of this formalism). Within this scheme, we 
find that even two-electron states can be confined to the dot. Figure 3 shows the ground 
state density for an interacting two-electron dot characterized by a radius of 15 units and 
the same spin-orbit parameters of Fig. 2. The density displays a circular ring-like shape 
although the spin-orbit term breaks the rotational invariance. As an additional check of 
this bound two-electron state, we have also solved the Hartree-Fock equations, finding a 
similar total energy and density distributions also confined to the dot region. Extending the 
calculations to treat three or more interacting electrons with the same spin-orbit term does 
not lead to confined states. 

In summary, we have proposed a mechanism to create electronic spatial confinement 
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based on the spin properties of III-V semiconductor structures and, more specifically, on 
the Bychkov-Rashba spin-orbit interaction and its strong dependence on the barrier layer 
composition. The eigenstates confined to these dots show common characteristics with 
conventional semiconductor quantum dots such as a discrete level structure of the single- 
particle states and the possibihty of confined two-electron states. The peculiar origin of 
the confinement reflects in the spin properties of the single-particle eigenstates. Not any 
material is suitable to create this type of structures since, as stated above, Rashba spin- 
orbit intensities of the order 10~^ eVcm and a range of tunability of the same magnitude 
are required. In this sense, InGaAs-based wells [7] could be enough to produce the confining 
effect at mesoscopic scale. However, since the bulk-intensity value Ag has a strong relevance 
in the determination of the minimum conditions for confinement, InAs-based wells reveal as 
the best candidates to fabricate these dots. This is because of their large built-in Rashba 
intensities, as large as 2 x 10~^ eVcm, and their tunability range of the same order [8]. 
This work was supported by Grant No. BFM2002-03241 from DGI (Spain). 
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FIGURES 
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FIG. 1. Structure of the bound states as a function of the generalized angular momentum j 
for a quantum dot having Ag = 0.2 H*aQ ~ 2 x 10~^ eVcm, Aj = OAH*aQ ~ 4 x 10~^ eVcm and 
a nominal dot radius Rq = 10 effective atomic units (lower scale). Solid lines represent the bulk 
bands corresponding to the surrounding material characterized by Ag (upper scale). 
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FIG. 2. Upper panels: radial spinorial components corresponding to the lowest energy pair 
of conjugated states. Middle panel: two-dimensional representation of the probability density 
corresponding to the above Kramers conjugated states. Lower panels: radial profile of the spin 
density distribution corresponding to the same conjugated states. All coordinates are expressed in 
terms of atomic effective units of distance 
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FIG. 3. Two-dimensional representation of the LSDA density corresponding to an interacting 

two-electron dot characterized by Ag = 0.2 H*aQ ~ 2 x 10~^ eVcm, Aj = 0.4 H*a^ ~ 4 x 10~^ eVcm 
and a dot radius of 15 units. Atomic effective units of distance are used to represent the Cartesian 
coordinates. 
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